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We present a new approach to the problem of Bloch electrons in magnetic ( some- 
times called Azbel-Hofstadter problem) field, by making explicit a natural relation 
between the group of magnetic translations and the quantum group Uq{sl2). The 
approach allows us to express the "mid" band spectrum of the model and the Bloch 
wave function as solutions of the Bethe-Ansatz equations typical for completely inte- 
grable quantum systems. The zero mode wave functions are found explicitly in terms 
of g-deformed classical orthogonal polynomials. In this paper we present solution for 
the isotropic problem. We also present a class of solvable quasiperiodic equations 
related to Uq{sl2). 



I. INTRODUCTION 

The peculiar problem of Bloch electrons in magnetic field 0, 0, [|], often emerges 
in various branches physics. Every time presenting a new face to describe another physical 
application. The Hamiltonian of a particle on a two dimensional square lattice in magnetic 
field is 

<ra,m> 



II e^^-™ = e^* (2) 

plaquette 

where $ = 27r^ is a flux per plaquette, P and Q are mutual prime integers and tn,m is a 
hoping amplitude between the nearest neighbors. In the most conventional Landau gauge 
Ax = Af^^fi+i^ = 0,Ay = ^rix the Bloch wave function is 

7A(n)=e^^>„^(fc), ^n = i^n+Q (3) 
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where Ux = n = 1, ...,Q is a coordinate in the magnetic celL With these substitution the 
Schrodinger equation turns into a famous one- dimensional quasiperiodic difference equation 
("Harper's" equation): 

tx{e'''^i^n+i + e-'^^i)n-i) + 2tj, cos{ky + = E^j^ (4) 

The spectrum of this equation has Q bands and feels the difference between rational and 
irrational numbers 0. The beuaty and complexity of this problem is the hierarchy of the 
spectrum. If the flux is irrational, the spectrum is singular continuum - uncountable but 
measure zero set of points (Cantor set) 0]. Themultifractal behaviour ( see for a review) 
of the spectrum is a long standing theoretical chalenge. 

The Harper's equation (i.e the Shrodinger equation for the Bloch particle in the Landau 
gauge) formally describes one dimensional electrons in quasiperiodic potential and also one 
dimensional quasicristal. Depending of the stength of the potential t^/ty (it corresponds to 
anisotropy in hopping amplitude of a particle in magnetic field) the Harper equation de- 
scribes localization- delocalization transition in one dimensional incommensurate potential. 

The list of others applications may be continued. 

Recently it has been conjectured that the symmetry of magnetic group may appear 
dynamically in strongly correlated electronic systems [|ll], [12 . 



In this paper we made explicit a long time anticipated connection of the noncommu- 
tatative geometry in magnetic field and the structure of the Quantum Group Uq{sl2) and 
therefore with Quantum Integrable Systems |T^. We present the Bethe- Ansatz algebraic 



equations for the spectrum of isotropic problem. Although we do not solve the Bethe- Ansatz 
equations here, we hope that they provide a basis for analytical studying the fractal prop- 
erties of the spectrum. Although, the Bethe Ansatz solution is available for an anisotropic 
case, in this paper we present only isotropic solution tx = ty = 1. Anisotropic case wil be 
considered elsewhere. 

The spectrum of our problem is complex but not chaotic. Instead it is governed by 
the Quantum integrability. It may happened that the Bethe Ansatz technique will be a 
powerfool tool to study general properties of multifractality. In fact multifractality has been 
already observed in quantum integrable models [|10|- However, this aspect of integrability 
has never been seriously developed. 

Let us stress that we apply the Bethe Ansatz technique to Quantum Mechanics - not to 
a Quantum Field Theory. Therefore we found convinient to use the theory of representation 
of the Quantum Group, instead of traditional Bethe Ansatz techique. 

The Harper's equation is a member ( perhaps the most interesting but at least the most 
famous) of class of solvable quasiperiodic equations related to Uq{sl2). It will be briefly 
described below ( for a more extended discussion of algebraization of quasiperiodic equations 
see Ref. [0). 

To ease the reference we state the main result now: 

It is known that due to the gauge invariance, the energy depends on a single parameter 
A = cos{Qkx) + cos{Qky) 0. We find that the spectrum at A = ( "mid" band spectrum 
) is given by the sum of roots zi 

Q-i 

E = iqQ{q-q-^)Y^Zu (5) 
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of the Bethe-Ansatz equations for the Quantum Group Uq{sl2) with 



n '=1 (6) 



q = e-^", (7) 

Another version of the Bethe-Ansatz equations is presented in Sect. 7. An ambitious problem 
of extending the Bethe Ansatz solution forarbitrary momenta k^, ky is beyond of the scope 
of this paper. The paper is organized as follows. In Sect. 2 we present the Schrodinger 
equation in two other gauges in which the quantum group structure of the model is more 
transparent. In Sect. 3 we review some necessary facts about the quantum group Uq{sl2) 
and its representations. The relation between the quantum group and magnetic translations 
(^) is revealed in Sect. 4 The Hamiltonian ([^) can be rewritten entirely via the quantum 
group generators . We give two equivalent ( dual) representations of the Hamiltonian as a 
linear or quadratic form in the quantum group generators. In Supplement A we consider 
a general quadratic forms in Uq{sl2) generators and introduce an "integrable " class of the 
second order "quasiperiodic" equations related to the quantum group. In Supplement B we 
show that a general quadratic form in quantum group generators is atrace of monodromy 
matrix of some integrable models with nonperiodic boundary conditions. In Sect. 5 we apply 
the functional Bethe Ansatz to obtain the "mid" band spectrum of the model. A promising 
connection between the zero mode wave functions and g-deformations of certain classical 
orthogonal polynomials is discussed in Sect. 6. 



II. MAGNETIC TRANSLATIONS AND GAUGES INVARIANCE 

The wave functions of a particle in a magnetic field form a representation of the group 
of magnetic translations : let generators of the translations be 

T^(i) = e^^M+M \ i><i + pi\ (8) 

They form the algebra 

rp^ rji — \ T^^T^ ^ — nxmrp ^ 

fl J- —jli J-n-'-m, q -'-n+mi 

TyT^ = q'^gTy, TyT_^ = q~'^T_^Ty (9) 

with q = expin^. The Hamiltonian (|I|) therefore is 

H = T, + T_, + Ty + T^y (10) 

Integrability of the problem can be seen from the first glance. It has been observed (see 
e.g. 0) that the spectrum depends only on one parameter A. That means that there is a 
parametric family of Hamiltonians 

H{u) = uT^ + M"^T_^ + vTy + v'^T^y (11) 
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with the same spectra at 

yQ^ik^ u-Q^-ik^ + w'^e*^^ + f ^'^e"*^^ = const 

. The mid band spectrum A=0 is a very spetial point where the Bethe Ansatz is especially 
simple. We shall consider only this point below. The useful way to explore this idea is 
representation teory of quantum groups. 

We found two specific gauges for which the quantum group structure is more transparent. 
In these gauges (we call them "regular") the wave function has the form 

Q-l 

= n i^O - Zm) (12) 

m=l 

where qo = q or qo = q^ and are the roots of the Eq.(8) and are independent on n. The 
Landau gauge of the Introduction is not "regular" in this sense and is not convenient for 
revealing the quantum group structure of the model. 

1. Modified Landau gauge 

Consider the gauge 

= -Ay = -(^n^ (13) 
which we refer as a modified Landau gauge. Now he Bloch wave function is 



Ak'n 



<PnAk'), <Pn = (pn+Q (14) 



where k' is the wave vector (different from that in (^)). The Schrodinger equation turns in 
to 

e'^'^-'^^cPn+i + e-*^i+**("-iVn-i + 2 cos(A:; + n$)0„ = E0„ (15) 

(we drop the argument in 0„). 

The "mid" band spectrum (A = 0) corresponds to the values 

^=(^<l>(g-l),7r) (16) 
in the Landau gauge and is given by 

k' = iO,n) (17) 

in the Modified Landau gauge. 

At the "mid" band all the values of k' on the line expi{k'^ ~ ^'y) ~ ~^ physically 
equivalent i.e. connected by a gauge transformation. We have chosen expz(fc^ + k'y) = —1 
indO). 

2.ChiraI gauge 

Consider the chiral gauge defined by: 



in wich the Bloch wave function takes the form 



Xin) = e^^Xn{p) (19) 

where p = {px,Py),n = Ux + ny,p± = {p^ -^Py)/2 are hght cone coordinates and momenta. 

This wave function is defined in two magnetic cells because n runs now from 1 to 2Q. 
Accordingly, the light cone momenta are confined to the half of the Brillouin zone [0, vr/Q]. 
The equivalent form of the Harper's equation is 

2e-i*-P+ cos(^$n - - P-)Xn-i = Exn (20) 

The wave function Xn is 2(5-periodic. This doubling of period in comparison with (^) is, 
of course, artificial. Although the coefficients in (pO]) have period 2Q there exists a simple 
transformation of Xn which makes them Q-periodic: 

Xn = exp(^n(n - 2))e„ (21) 
This new wave function ^„ is Q-periodic and satisfies the equation 

(e3*-^P- + e-'*"+T*-*Pj/)^„_i = (22) 

The "mid" band spectrum corresponds to p = (|7r, ^tt + in the chiral gauge (compare 
with (p^ ). Again all the values p = {\tt,S) for any real 6 are physically equivalent. We 
choose p at the "mid" band to be 

(23) 

The relation between the three gauges are given in the Appendix A. 



III. QUANTUM GROUP 



The algebra Uq{sl2) ( a q-deformation of the universal enveloping of the s/2) is generated 
by the elements A, B,C, D, with the commutation relations |2^, 



AB = qBA, BD = qDB, 
DC = qCD, CA = qAC, 

A2 _ r)2 

AD = 1, [B,C] 



q-q 



-1 



(24) 



We shall take the deformation parameter g as a root of ±1 of degree Q: q = exp{i7r^) where 
P and Q are mutually prime integers. The central element of this algebra (for arbitrary q) 
is a g— analog of the Casimir operator 
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When g is a root of unity some additional central elements appear. 

In the classical limit q — *• 1 + |$, the quantum group turns to the 5^2 algebra: 
{A - D)/{q - g-i) ^ S3, 5 ^ S+, C ^ S.,c^ + 1/4. 

The commutation relations (p^ are simply another way to write the interwining relation 
for the L - operator: 

R{u/v){L{u) ® 1)(1 ® L{v)) = (1 ® L{v)){L{u) ® l)R{u/v) (26) 

with the trigonometric i?— matrix 



R{u) = ]^{q + l){u-q \ ^) + i(g-l)(?i + g \ ^)a^®G^ + {q-q ^){a+ 



o-_ + o-_ (X) (J I 



(27) 



satisfying the Yang-Baxter relation (a^ are Pauli matrices; a± = (ai ± ia2)/2). Generators 
A, B, C, D are matrix elements of the L- operator 



L(u) 



r ukA-u-^k~^D 

B 



C 

q-q-^ 



(28) 



Here u is the spectral parameter and k is an additional parameter (rapidity at the site). 
Note that the R- matrix is the L- operator in the spin 1/2 - representation. It is given by 
the same matrix (^Sf ) for k = g^/^ with elements: A = gz'^s^ ]j _ _ ^-^^ (j _ 

Irreducible finite 2j + 1 dimensional representations can be expressed in the weight basis 
where A and D are diagonal matrices: A = diag [q^ , q"-'). An integer or half- integer j is 
the spin of 

the representation. The value of the Casimir operator (^) in this representation is given 
by the q- analog of (j + 1/2)^. 
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-0 = ) = b + 1/2]^ (29) 

The representation can be realized by difference operators acting in the space of polyno- 
mials "^{z) of degree 2j . 

A^{z) = q-^^{qz), D$(;z) = q^^^^q-^z), 
5^(z) = z{q - q-^)-^ (g2%(g-i^) _ g~2%(g^)) 

C^{z) = -z'\q - q-^)-^ {^{q'^z) - ^{qz)) (30) 

Then ^q{z) = 1 is the lowest weight vector whereas ^2j{.z) = z"^^ is the highest weight 
vector,i.e. C^!o{z) = 0, B'^2j{z) = 

We call this series of representations "regular". They are a smooth deformation of the 
representation of the s/2 algebra by differential operators: 
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In addition, in the special dimension Q there is three parametric family of representations 
having, in general, no lowest and no highest weight. Sometimes they are called cyclic or 
unrestricted representations. These representations have no classical finite dimensional limit. 
They can be written in the Weyl basis p9|, |2^, Let X and Y satisfy 

qXY = YX (32) 

Then the representation of Uq{sl2) is 

A = X, D = X-\ 
B = {bX + bX-^)Y-\ 

C = Y{cX + cX-^), (33) 

and is characterized by parameters b,c,b,c obeying the conditions qbc = q'^bc = —{q — 
g,-i-j-2 rjj_^g value of the Casimir operator (|25| ) depends on the parameters and is equal 
to cb + cb — 2(g — g~^)~^. Comparing it with (pQ]) we find that at q%lb = c/ c = ^q (for 
P-odd, even) cyclic representations belong to the "regular" series with g^-'^^ = =pl for P - 
odd (even) and the Casimir operator (^9]) is 

Co = -4(g - g"^)"^ forP- odd 

Co = 0, for P — even (34) 



IV. MAGNETIC TRANSLATIONS AS A SPECIAL REPRESENTATION OF THE 

QUANTUM GROUP 

For a given value of the Bloch wave vector k the dimension of physical Hilbert space of 
our problem is Q. A representation of dimension 2j + 1 = Q of the quantum group Uq{sl2) 
for q = exp{i7i^) naturally acts in the space of Bloch states in the magnetic field. Therefore 
magnetic translations and the Hamiltonian (1),(5) can be expressed through generators of 
the quantum group. Let us set 

Ty = e'P^YX-^, 



YX, (35) 



Then, the Hamiltonian (|1^) could be expressed as a linear form in quantum group gen- 
erators. At the middle of the band (an integrable point) 

^i{p.-Py) = (36) 

we may choose 

b = q\q-q-')-\ 
c={l-q')-\ 

c = e5(P«-P-)(l-g2)-i, (37) 
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Then B and C given by (^) form a representation of the quantum group with the value 
of the Casimir operator Eq.(|3^) wich corresponds to the regular representation i.e. with 
the highest and the lowest weight. Using (|33| ) and (|35| ) we may identify the quantum group 
generators and magnetic translations 

T_, + T.y = {q- q-')e-^^^Py+P^'>B, 

T_,Ty = ±qD^ (38) 



The Hamiltonian now acquires the form 

H={q- q-^){^e^^Py+P-^C + e-5(P''+P-)B) (39) 

In fact at the integrable point (^) the physical states and their energies ( midband states) 
do not depend on the value of exp ^{py + p^). It simply can be gauged away. As a result for 
the midband states in the isotropic case the relation between magnetic translations and the 
quantum group is as follows 

H = i{q - q~^){C ± B) (40) 
The Hamiltonian (EO) can be written as the trace of a modified L- operator (ESF) L[u). 



Say at P odd 

H = TrL{u) = TrL{u)ai (41) 
Modified L-operator also obeys the intertwining relation (EBI) S 



Realization of the quantum group in terms of magnetic translations is not unique. There 
are different realizations wich may be used for different physical applications. Below we 
present, another "dual" realization. Say for an odd P we have: 

T_., + T_y = -i{q-q-^)q--^BD, 
T, + Ty = -i{q - q-')q--^CA, 

T_yT^ = q-^A\ T_,Ty = qD^ (42) 



that also forms a representation of Uq{sl2) with the same value of the Casimir operator (^91) 
and with the same dimension Q. Now the Hamiltonian turns into quadratic form in terms 
ofthe Uq{sl2) generators: 

H = -i{q-q-^)q--2{CA + BD) (43) 

Two different forms of the Hamiltonian are in fact gauge equivalent: they correspond to 
the two choices of gauge discussed in Sect. 2. To stress ambiguity of representation of the 



1* 



L-operator of such kind is used in the sine-Gordon model |22]. We are grateful to L.Faddeev 



and A.Volkov for pointing our attention to the Ref. ]22] 
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quantum group by magnetic translations we notice that the Hamiltonian of our problem can 
be also represented as a quadratic form of another quantum group Uq2[sl2) 

H = CA- BD + q^BA - q^CD (44) 



Diagonalization of this form leads to the same Bethe Ansatz equation as the form 

In the Appendix B we show that the Hamiltonian ( ^31) represented as a quadratic form of 
the quantum group generators can also be understood as a trace of a quantum monodromy 
matrix, though of a more complicated nature - the monodromy matrix for a system with 
non-periodic boundary conditions. 

Since the main ingredients of the quantum inverse scattering method have appeared in 
the problem there is a strong evidence that the Hofstadter Hamiltonian can be actually 
diagonalized for arbitrary Q. i.e. its energy spectrum may be expressed by solutions of 
the Bethe equations. The ambitious goal of solving the problem for arbitrary parameters 
Px.Py and tx/ty is beyond of the scope of this paper. Instead we shall present a solution 
for the fixed values of the parameters when representation of the quantum group (expressed 
through magnetic translations) belong to the "regular" series and, therefore, has the lowest 
and the highest weight. This point corresponds to the middle of band and isotropic hoping. 
In this case the so-called functional Bethe Ansatz is applicable. 



V. FUNCTIONAL BETHE ANSATZ 

Tthe Hamiltonian of the Bloch particle in a magnetic field in the midband point of the 
spectrum A = can be expressed as a polynomial form in the quantum group operators in a 
regular representation. Let us use this advantage and replace generators by their functional 
realizations. 

Let us first consider the representation in which the Hamiltonian is a linear form(^D|). 
Substituting the functional realization (|30| ) in the Hamiltonian ( ^OD we obtain a difference 
equation for a meromorphic function '^{z) 



+ qz)^{qz) - i{z-^ + q-^z)'^{q-^z) = ^^(z) (45) 



This equation can be easily obtained directly from the original Harper equation in the chiral 
gauge (|20| ) without a reference to the quantum group. Indeed, Eq.(^5D is an extension of 
the (^) at p = (f , f ) to the whole complex plane provided 

Xn = ^(g") (46) 

However the miracle which would be very hard to see without the quantum group is that 
the energy in the extended equation (|45| ) does not depends on z. 

Now it is easy to get a family of Hamiltonians which have the same spectrum. Let us 
set z = q^u and '^{q^u) = Xn{u). Then, all Hamiltonians 

H{u) = i{q-^u-' + q''+'u)Xn+i{u) - t{q-''u-' + q''''u)Xn-i{u) = Exn{u) 

have the same spectrum and correspond to A = 0. The original Harper's Hamiltonian (^ ) 
is Hil). 
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Another miracle is that we know in advance that the function "^{z) is a polynomial of 
degree Q — 1: 

^iz)=Y[iz~z^) (47) 

m=l 

Note that last fact is not valid in an arbitrary gauge, in particular in the Landau gauge. 

The similar arguments work in case of the modified Landau gauge. The functional 
realization of the quantum group Hamiltonian ( ^31) gives the difference equation 

z-^^{q^z) + q-^z^{q-^z) - {z + z-^)^(z) = E-^{z) (48) 

with a polynomial solution. It is an extension of the Harper equation in the modified Landau 
gauge for k' = (0, vr) can be extended to the whole complex plane with 

0n = ^{q'n (49) 



The most suitable method to solve the spectral problems (^ , ^81) or is the so called 
functional Bethe-Ansatz. 

Let us plug ( ^ in (0,|5) and divide both sides by "^{z). 

r{z-' + qz) n ff^^ 

-^[z-' + q-'z) \{ ^^^^^ = E (50) 

m=l,mj^l ™' 

The l.h.s. of this equation is a meromorphic function, whereas the r.h.s. is a constant. To 
make them equal we must null all residues of the l.h.s.. They appear at z = 0,a.t z = oo 
and at z = Zm- The residue at z = vanishes automatically. 

The residue at 2; = oo is —iq^ + iq"^. Its null determines the deqree of the polynom. 

Comparing the coefficients of z^~^ in the both sides of Eq.(^),we obtain the energy 
given by Eq.(^ advertised in the Introduction. 

Finally, annihilation of poles at z = Zm gives the Bethe-Ansatz equations (^) for roots of 
the polynomial (0). Here we write them in a more conventional form. Let zi = exp {2ipi), 
then 

cosh(2(^; - if) ^ ^ '^-i sinhjipi -ipm + ij) . . 

COsh(2y?; + if) ^=i^rn^i sinh(^; -^m- ^f ) 

The Bethe equations for the quadratic realization of the Hamiltonian (^3]), can be obtained 
via similar method. They are 

zf = q^ n / = l,...,g-l (52) 

m=l,m^l 1 ^rn 

The energy spectrum is again proportional to the sum of roots: 

Q-i 

E = -q(q - g-i) ^ zi. (53) 
1=1 

Inspite of the apparent difference eqs. (^) and (|) must be equivalent to the eqs. ( |52D and 
(0). The former is useful tro describe the middle of the spectrum, whether the latter is 
good for the bottom of the spectrum. 
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VI. Q-ANALOG OF ORTHOGONAL POLYNOMIALS AS EXACT ZERO MODE 

WAVE FUNCTIONS 



There is an intriguing connection between the wave function of the zero energy state 
and q-generahzation of the classical orthogonal polynomials The q-analoq of orthog- 
onal polynomials (so called Askey- Wilson polynomials) are the most general polynomials 
orthogonal on a descrete support. They depends on four parameters (except q) and satisfy 
a q-analog of the differential hypergeometric equation. It is the difference equation 

A{z)Pn{q'z) + A{z-')P4q''z) - iA{z) + A{z''))P4z) = {q-'^ - 1)(1 - abcdq'"^-') Pn{z) (54) 



where 

_ (1 ~ ~ M(l - c^)(l - dz) /c-c-N 

(1 — ^ )(1 — q^z'^) 

and a, b, c, d are parameters and n is the degree of the polynomial. Note, that P„ is a Laurent 
polynomial in z and usual polynomial in 2; + z~^ of degree n. Recently, the Askey- Wilson 
polynomials appeared in connection with representation of the quantum group. It has been 
found that they are closely related to 6 — j -symbol of Uqi^sh) f^^. They also appeared in 
"g-harmonic analysis" to be spherical functions on the quantum SL(2) group Some of 
them are closely related to our problem. 

Choosing c = —d = q, a = b = we arrive at the equation for the continuous g-Hermite 
polynomials ||3^ H^^ 



H^\q'z)-z'H^\q-'z) = q-^-{l-z^)H^^ 



(They are called continuous because of their orthogonality on the unit circle with a contin- 
uous measure). 

It is clear from (|48|) that for an odd Q g-Hermite polynomial of the order n = (Q — l)/2 
yields zero energy solution of our problem: 

*(^=°)(z.) = .(«-i)/^/7g_,)/,(.). (56) 

The explicit form of the g-Hermite polynomials is 

where the standard notation 

n-l 

(«;?)n = ]]_{l-aq^) 

1=0 

is used. 

Another choice is c = —d = q, a = —b = q and then the replacement of g by g^/^ gives 
the q-Legendre equation 

^-^^Pi'Kl^) + J^Pi'\l''^) = + l''^')Pi'\^) (58) 

A. A. 
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Comparing with the Eg. (^5]) we conclude that the zero mode solution is given by the con- 
tinuous g-Legendre polynomial 



Their explicit form is 



p(l)(^z''j = ?)n(g'^'^^; Q)m{<l^^'^] ll)n-m ^2m-n /gg\ 

171=0 

n—m 



VII. CONCLUSION 

We have showed that motion of a particle in a periodic potential in a magnetic field 
has posses a structure of the quantum group. We found the Bethe Ansatz equation for 
the midband spectrum and the wave function for a particle on a square lattice. We also 
presented a class of quasiperiodic equations related to Uq{sl2) which can be solved in a 
similar way. The most interesting feature of that type of equation is the multifractality of 
their spectrum at P, Q — > oo, i.e when flux $/27r is irrational. The spectrum is complex but 
not chaotic. Quite opposite, it is determine by the quantum integrability. 

As we already mentioned the Bethe Ansatz solution is available (but not presented here) 
for an arbitrary strength of the potential of the Harper equation. It provides the basis to 
study Anderson localization-delocalization transition of electrons in quasiperiodic potential. 

As the best of our knowledge even elementary questions regarding the Hofstadter problem 
at irrational flux remained unclear. For example the low temperature thermodinamics, dis- 
persion of excitations, conductivity etc. are not known. In all previous examples of quantum 
integrable models all these quantaties have been found by solving Bethe Ansatz equations 
at large Q. In this limit dispersion of bands becomes negledgible. Therefore, the lack of 
solution away from the midband is therefore is not an obstacle for study multifractality and 
other physical properties. The Bethe Ansatz solution does not give an advantage at finite P 
and Q. Just contrary, up to Q ~ 6000 direct diagonalization of the Hofstadter Hamiltonian 
is more effective. However, as usual the Bethe Ansatz solution becomes a powerful tool at 
Q ^ oo. Perhaps the most interesting task is to solve the Bethe Ansatz equations in this 
limit. The strategy of solving the Bethe equations is well known: at large Q the roots zi 
form dense groups ("strings") and can be described by their distributions. The algebraic 
equations then replaced by the system of integral equations for the distribution functions of 
strings. This program is in progress. 

Another ambitious problem is to generalize the Bethe ansatz solution to the whole Hilbert 
space of the model. 



VIII. SUPPLEMENT A:A CLASS OF DIFFERENCE AND DISCRETE 
EQUATIONS RELATED TO THE QUANTUM GROUP 

A genral quadratic form in generators of the quantum group provide a class of difference 
and discrete operators of the second order which allows complete or partial algebraization. 
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A problem of the Bloch particle in magnetic field is a particular case (the most physically 
valuable, perhaps) of the general class. 

The idea of algebraization is simple. Consider a regular representation of the quantum 
group with the lowest and the highest weight. The space of this representation can be can 
be realized by the ring of polynomials of degree 2j. Any form in generators A, B, C, D 
preserves the space of polynomials. Therefore eigenfunctions of a form, which by means 
of the representation ( pil| ) is a difference operator, are polynomials as well. Vice versa is 
probably also true - one may look on the quantum group as an algebra of difference operators 
which preserve the ring of polynomials. 

Quadratic forms of the "classical" s/2 as a class of second-order diffferential equations 
solvable in polynomials (so called "quasi exactly solvable" models of quantum mechanics ) 
have been extensively studied I^G] (see also |^ for a revew) in the last few years. They 



include second order differential equations for all classical orthogonal polynomials and a set 
of Schrodinger operators with solvable potentials. Attempt of q-generalization for difference 
functional equations (in fact very closed to our approach) has been made in Ref. i. 



Consider a general quadratic form (in the Appendix B we show that it can be consider 
as a trace of monodromy matrix of some integrable model) 

G = aA^ + dD^ + {q- q-^)ic2CA + b2BD + hBA + c^^CD) + (g - q-^fih^B^ + CiC^) (60) 

with a set of parameters a, d, q, hi {i = 1, 2, 3). 

Under the representation ( pO| ) the quadratic form is a difference operator. Let us consider 
its spectral problem. 



where 



viz] 



G^{z) = a{z)^{q'z) + d{z)^{q-'z) - v{z)^{z) = E^{z) (61) 

a{z) = b,q-^'^'z^ - hq-'^'z + aq-""' + Csg^^z"' + Ciq-^"^ (62) 

d{z) = hiq^^'^z^ + h2q^^z + dq^^ - c^q^ z-^ + Ciqz'^ (63) 

= {q + q-^){hz^ + ciz-^) + {C2q-^ - c^q^)z-^ + (foag"'' - hq')z (64) 



If all 61^2,3 = the difference equation (|6ll) is known as g-hypergeometric equation (see e.g. 
p^). Their solutions are certain g-Jacobi polynomials - a degenerate case of Askey Wilson 



polynomials (M) [pq] . At the "classical" limit g — 1 they reproduce classical orthogonal 
polynomials. The difference equation ( pT| ) in it general form in the classical limit leads to 
Mathieu and Lame equations. One may call the equation (|6T|) g- analog of the Mathieu 
(Lame) equations. 

From the theory of representation of the quantum group, we know that this equation 
has (2j + 1) polynomial solutions of degree no greater than 2j. 



^* We are indebted to A.Turbiner who informed us about this work 
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N 

^{Z) = n - ^rn) (65) 
m=l 

The class of "integrable" equations can be extended if we apply a "gauge" transformation 

^{z) ^ f{z)^{z), a{z) ^ a{z)f{z)/f{q'z), d{z) ^ d{z) f (z) / f (q-^) , v{z) ^ v{z) 

A class of "integrable" discrete equations can be obtained from (^) by setting z = 

a„^„+i + - = n = 1, Q (66) 

Solutions of the discrete equations are given by discrete polynomials 

= n - ^rn) (67) 

m=l 

with the same roots. 

The functional Bethe-Ansatz described above determines algebraic equations for roots 
of the solution (|65|) of eqs. (^) or (|66| ) i. We already applied this method for the Harper 
Equation. For complete reference we repeat it for a general case below. First of all, let us 
plug (^) in (|6lD and divide both sides by "^{z). We get 

N 2 N -2 

n '-^^ + d{z) n - v{z) = E (68) 

Z Z^ Z Zra 

If at least one of the coefficients ci, C2, C3 is nonzero, there are two different cases: 

(i) at least one of 61, 62, &3 is nonzero, 

(ii) All y s are zero -the quadratic form (|60|) includes only A, D and C (generators 
of the Borel subalgebra of Ug{sl2)). That type of difference equation (|61|) is known as 
q-hypergeometric equation. Their solutions are certain g-Jacobi polynomials - a class of 
polynomials orthogonal with a discrete measure (see e.g. pl[). 

Let us consider the case (i) first. The l.h.s. of (|68D is a meromorphic function, whereas 
the r.h.s. is a constant. To make them equal we must cancel all the singularities of the l.h.s. 
They appear at singular points of a{z), d{z) and v{z) (double and simple poles at z = and 
z = 00) and at z = Zm- 

The singular part at 2 = vanishes automatically. 

Vanishing of the singular part at z = 00 

^^(^2JV-4,+i ^ ^-2^+4,-1 _ ^ _ ^-ly ^ b2{q-'''^'' - q-')z + h{q^ - q^''-^^)z + 

N 

+z{q - q-'Mq'"^-'^ - g-^^^^^) ^ z^ (69) 

m=l 



3** 



We note that this method is most close to the approach of algebraization of quantum me- 



chanical problem suggested by Ushveridze |26] 



14 



determines the degree of the polynomiah = 2j. 

Comparing the constant terms in the both sides of (|6^) we find the energy spectrum: 



N N 
n<m m=l 

(70) 

Finally, annihilation of poles dX z = Zm gives the following Bethe-Ansatz equations 

4^ = .^^- 11 / = 1,...,2,. (71) 

The Bethe equations is a system of 2j algebraic equations. It must have exactly 2j + 1 
solutions corresponding to different eigenfunctions. In the case (i) all of them are polynomials 
of one and the same degree 2j. 

Similar arguments are applicable in the case (ii) . The difference is that the l.h.s. of (|68D 
is now regular at ^ = 00 from the very beginning and the condition of vanishing of ( |OUD does 
not bring a restriction on degree of the polynomials. The Bethe equations are valid for any 
iV < 2 j + 1 

^ = q-^ n l = l,...,N. (72) 

Apparently, for each degree there is exactly one such polynomial. This means that for each 
A^ < 2j + 1 the Bethe equations ([72|) must have exactly one solution. This fact is far from 



obvious when we look at (|7^). Anyway, we derived equations (|7^) for zeros of g-orthogonal 
polynomials. It seems to be interesting to review the theory of orthogonal polynomials from 
this point of view. 

The main property of the g-hypergeometric equation (the case (ii)) that makes it "solv- 
able" is triangularity of the matrix connecting the original basis 2;™ (m = 0, 1, 2j) with 
the basis formed by eigenfunctions of the operator G. In particular, this leads to a very 
simple structure of the spectrum of G: 

E^v = ag2^-2^^ + rfg2^-2^ Ar = 0,...,2j. (73) 



IX. SUPPLEMENT B:QUANTUM INTEGRABLE MODELS WITH 
NON-PERIODIC BOUNDARY CONDITIONS AND HOFSTADTER PROBLEM 

A general quadratic form in quantum group generators is related to a quantum magnetic 
chain with one site and non-periodic boundary conditions. 

To begin with, let us give a brief summary of the formalism treating the integrable 
systems with boundaries. The boundary conditions of an integrable model are determined 
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by c- number 2x2 matrices K^{u) and K^{u) depending on spectral parameter and satisfying 
the "reflection equations" |l27|l . 



R{u/v){K_{u) ® l)R{uvq-^){l K_{v)) = (1 ® K_{v))R{uvq-^){K_{u) ® l)R{u/v) 
R{v/u){KX{u) ® l)R{{uvq)-^){l ® KX{v)) = (1 ® KX{v))R{{uvq)~^){Kl{v) ® l)R{v/u) (74) 

{t means the transposition) with the -R-matrix p7|. Each solution of the "reflection equations" 
specify a boundary condition consistent with integrability. Solutions for and for K_ are 
related 

K+{u) = K'_{u-^) (75) 

The monodromy matrix for an integrable model with non-periodic boundary conditions is a 
quadratic form in the monodromy matrix L{u) of a model with periodic boundary conditions 

m 



T{u) = L{u)K^{u)a2L\u'^)(T2 (76) 

The trace of the monodromy matrix 

t{u) = Tr{K4u)T{u)) (77) 

forms a commutative family [t{u),t{v)] = 0. 

It is known that for the models with the trigonometric i?-matrix there is a 3-parametric 



family of boundary /^-matrices [pQ], |31 



K(u) = (^^^ ^u-qsu 1) f3{q ^u'^ - qu 2) \ 

_ ( \{qtu - q-H-^u-^) ^{qu^ -q~^u-^)\ . . 

+ iy{qu^ - q-^u-^) -\{t~^u - tu-^) ) 

where a, /3, 7, s and A, /x, z/, t are arbitrary parameters characterized the boundary condi- 
tions. Substituting the L-operator (^) into ( |75D we find matrix elements of the transfer 
matrix 

^ a{q-'s-'u- qsu-') , ^ _ _ ^ 
[q-q-^y 



q u — qu 



q-q- 



- — {-fku^^CD - -fk-^uCA - fSkuAB + (3k-^u-^DB) + a{su - s-^u-^)CB (80) 



q-'u'-qu~\Pik'A' + k^'D^-u'-u-') 

Ti2{u) = — ( l{q-q )C + 

q — q Q ^ Q 



+a{sk'^DC - s-^kAC)) (81) 



T21 M = I3{q-q )B 

q — q~^ q — q~^ 

+a{skBD - s-^k-^BA)) (82) 
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q u — qu 



{l3ku~^BA - (3k-\BD - ^kuDC + ^k-\-^AC) - a{q'h-\ - qsu-^)BC (83) 

q-q 

Summing diagonal elements (pO]) and (p3|) we find the trace of the transfer matrix t{u) 



t{u) = ^""^ ^ ""^^^ ^lly ^ + ^I3k'' - a\ts-^)A^ + {^^^ + z//3fc-2 - a\t-^s)D^ + 

{q - q'^){-{tk~^^l + q'h-^kav)CA + {q-H-^k-^\p + ska^ji)BD + 
(r^A;A7 + qsk-^av)CD - {qtkXp + s-^k-^afi)BA) - 

— (g — q^^Y(nPB'^ + i"jC'^) ] + c — numberterm (84) 

(here we replace the Casimir element (|2^) by it c-number value ) 

The trace of the transfer matrix r(-u) is a general quadratic form in A, B, C, D since it 
depends on 7 parameters - 3 in each i^-matrix and rapidity k. Indeed, the total number 
of coefficients of a general quadratic form is 10 but one of them is a common multiplier 
and another two contribute to the c — numberterm in ([5^) due to the two central elements 
{AD = 1 and the Casimir operator). 

As we showed previously the Hamiltonian of the Bloch particle in magnetic field is a 
particular quadratic form of the quantum group generators and therefore can be considered 
as an integrable model. 

APPENDIX A: A 

Gauges 1. The connection between Landau gauge and the modified Landau gauge ([I3D 
is very simple 

<j)n{k') = exp(yn(n - Q))Mk) (Al) 

k'x = kx + - ^, k'y = ky (A2) 

Both k'^ and k^ are defined modulo $ . 

2. The relationship between Landau gauge and the chiral gauge is a bit sofisticated 0We 
are indebted to A.Abanov who found this connection. They are connected by the Fourier 
transformation 

0™ = E e^*"'"^, (A3) 

n=0 



4* 
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where (pn is the wave function in the modified Landau gauge. The new function (/)„ obeys 
the equation 



He-^'y + e-**"+**+^'=^)0„_i = E0„ (A4) 
Comparing with (22) we can identify the momenta as follows: 
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Postscriptum: when this paper has been finished we have learned that L.Faddev and 
R.Kashaev succeded in extending the Bethe Ansatz solution for an arbitrary momentum 
and an arbitrary anisotropy. We indebt to L.Faddeev and R.Kashaev for informing us of 
their results. 
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